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•ABC TRACT 



An experimental study of the diffusion of helium in the wake of a 
circular cylinder was conducted in the GALCIT hypersonic wind tunnel 
at a Mach number of 5. 8. The cylinder was constructed of material 
having random porosity and was mounted with its axis perpendicular to 
the stream. The light gas was injected in small amounts and the thermal 
conductivity method was utilized to detect the concentration of helium in 
the air at points downstream. Problems in the utilization of the thermal 
conductivity method for low sample densities were overcome by suitable 
calibration. 

Flow in the wake of the cylinder was found to display character- 
istically similar behavior at a few diameters downstream, with respect 
to decay and spread of the concentration. Reynolds number similarity 
was established in the laminar case, but turbulent Reynolds number 
similarity may require reference to momentum thickness, which was not 
possible with the present data. 

Profile data was somewhat marred by a tunnel pressure perturba- 
tion, but many of the important conclusions were not affected. The pro- 
files appear to follow the theoretical Gaussian distribution in the similar 
region. 

The thermal conductivity method is quite promising as a means 
of tracing the diffusion of one binary gas constituent in another, as 
applied to hypersonic wind tunnel experiment. It will also serve in the 
analysis of transition and turbulence, and of the lateral spreading of 
the turbulent fluid into the rest of the wake region behind the bow shock. 
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I. INTRODUCTION 

The problem of the wake behind a blunt body is one of the oldest 
in the classical, incompressible aerodynamic literature. The practical 
aspects have centered around the effects of wake vorticity and circulation 
and the momentum defect determination of drag*#. A limited quantity of 
this classical literature has been directly concerned with the problems 
of mass diffusion in the wake. There exists, however, a great amount 
of theory and experimental data on jet mixing o£ incompressible flows 
and lower Mach number compressible flows. The jet mixing theory is 

of course directly analogous to wake theory in many respects, especially 

2 

in the subsonic case . Also, for the viscous-temperature effects of 
primary interest here, namely those of the diffusion of momentum, 
mass and energy, the jets and wakes are just free boundary cases of 
boundary layer flow for which there is extensive treatment in the 
literature^. 

A summary of the mixing theory and list of references to 1954 

Z 

has been contributed by Pai . This compilation indicates that the know- 
ledge of diffusion processes of incompressible flow is fairly complete, 
and consists for the most part in consideration of momentum-vorticity 
diffusion. The compressible subsonic theory is also well established, 
with the addition of some topics of heat energy diffusion. More recent 

contributions which are of particular interest in the present investigation 

4 5 

are the works of Townsend', and Schubauer and Tchen . Townsend 
proposes such ideas as that of the intermittency at outer boundaries of 
turbulent flows and discusses the wake behind a two-dimensional 



* Superscripts denote references at the end of the text. 
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cylinder at some length. Schubauer and Tchen collect and extend the 

ideas of Townsend and others in a well-integrated description of such 

flows, and include an excellent bibliography of the associated literature. 

An additional recent contribution to the specific field of interest here is 

6 

the book by Hinze . 

Literature on the supersonic diffusion processes in wakes cannot 
depend so heavily on the jet mixing literature, and very little has been 
done, even recently, in this more specific field. Also, the literature 
described indicates little concern with mass addition and chemical 
reactions, which have recently become highly important in tire hypersonic, 
high altitude flight regimes of advanced weapon design and space explora- 
tion. A number of theoretical and experimental papers have been written 
in recent years which have demonstrated clearly the importance of these 

physical- chemical considerations in the transfer of heat to bodies in 

7-15 

high Mach number flow media . A good summary of associated work 
7 12 13 

is given by Lees . Hubesin and McMahon , in particular, show that 
a light gas is most effective with respect to weight of material injected 
in reducing heat transfer rates. One concludes from such a literature 
survey that many hypersonic problems of practical nature will involve 
the injection of materials in light gaseous state into the boundary layer 
for cooling, which in turn may result in chemical reactions and nonuni- 
form mixtures of gases in the wake flow behind the body. Any complete 
study of the hype'r sonic wake must consider these phenomena. 

The wake flow behind a hypersonic body involves the usual con- 
siderations regarding afterbodies, structure, and control surfaces, but 
it has other implications as well which may become quite important. 

There is, for instance, the probable radar reflectivity and optical 
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emission properties of the material in the wake, which are of interest 
in the study of meteor wakes and in the tracking of reentry vehicles, in 
addition to other obvious military considerations. (Dr. Feldman and 
his associates at AVCO Research Laboratories are currently investiga- 
ting these problems. ) On the other hand, the study of the hypersonic wake 
is of current interest as a part of the present exhaustive investigation of 
the hypersonic flow field in general. 

The above considerations have led to the establishment at GALCIT 
of a program for systematic study of the hypersonic wake behind blunt 
bodies. Wake pressure measurements at a nominal Mach number of 
5. S have been undertaken by J. M. McCarthy, and temperature and hot- 
wire aaemometry measurements are being conducted by A. Bemetriades 

and C. F. Dewey. In a separate study, M. D. Coffin is continuing the 
13 

work of McMahon in the investigation of heat transfer with mass 

injection at the stagnation point of a blunt body. In his work. Coffin 

has developed an apparatus for the analysis of the concentration of one 

gas in another by the method of thermal conductivity. The method was 

16 

known to him through his contact with the work of Rush and For stall , 

1? 

and Forstaii and Shapiro , at the Massachusetts Institute of Technology. 

The present investigation was suggested by Captain Coffin (USAF) 

and Professor Lees, as a means of extending the GALCIT wake program 

by a diffusion study. The use of a tracer gas should help in understanding 

the chemical aspects of wake flow and diffusion in a binary mixture 

of gases. It was also envisioned that a light gas might be used as a 

tracer to determine turbulent processes in the wake. Although the 

1 8 

thermal conductivity method is a classical one in gas analysis , it 
has seen very little employment in the field of wind tunnel experiment. 
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In view of some of the limitations of other devices in the usual hypersonic 
wind tunnel apparatus for the measurement of the physical state of the 
gas, it has become evident that new methods for obtaining and correlating 
this type of data are required. A large portion of the discussion below 
will necessarily involve the setting down of simple procedures, 
corrections and considerations that have been found essential in the 
employment of the apparatus. It is hoped that inclusion of all relevant 
procedures will assist future investigators in avoiding the many small 
technical delays that interfere with a smoothly conducted experiment. 
Some improvements to the equipment which could not be incorporated 
in the present study will also be indicated. 

Thus, the present investigation is a study of the effectiveness 
of the thermal conductivity method for determining concentration of a 
tracer gas. At the same time it is an exploratory study of the utility 
of this tracer gas in uncovering the nature of hypersonic wakes behind 



blunt bodies. 
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II. EQUIPMENT AND PROCEDURE 
II. 1. General Experimental Method 

In accordance with the cited objectives the following experimental 
steps were planned: 

(1) Ejection of light gas from a porous body of simple geometric 
shape at a nominal Mach number of 5. 8 in carefully controlled small 
quantities. 

(Z) The extraction of samples of the wake at various stations in 
the tunnel by means of a small pitot tube. These samples were to be led 
to a thermal conductivity ceil which would compare concentration of the 
diffusing agent in the sample with a known reference gas. This procedure 
requires a calibration of the thermal conductivity of known mixtures of 
the two gases against cell readings. 

(3) The correction of the data for mean tunnel concentration, 
differences in pressure of the sample and reference cell, instrument 
errors, and any other effect that might alter consistency and reproduci- 
bility of the readings from test to test. 

(4) The comparison of wake pressure data taken for the ejection 
model against that of the solid model used for total head surveys, 
(McCarthy's model), to correlate results. 

(5) Experimental check of the two-dimensionality of the tunnel 
flow at the center of the model, and spanwise uniformity of the tracer 
gas ejection. 

(6) The reduction of the data to find what conditions of similarity 
may be found in the flow, and whether turbulent zones can in fact be 
defined by the method. 
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(7) The coir- pari son with available theory and related experiment. 

II. 2, Selection of Diffusion Gas 

The selection of a tracer gas for diffusion into the tunnel air 
that would best meet the requirements of the experiment was simple. 

If one checks the list of the common unobjectionable, relatively 
inexpensive laboratory gases, helium and hydrogen are the two most 
suited to the application and most sensitive to thermal conductivity measure- 
ment. Helium is 5. 97 and hydrogen 7. 15 times as conductive as air, but 

1 6 

hydrogen is inflammable . (See Table I. ) The choice of helium becomes 
obvious, and it is certainly one of the light materials that might be 
employed in surface cooling applications. The commercial helium 
utilized was 99. 9976 per cent pure, with traces of carbon dioxide, 
argon, hydrogen, nitrogen, and methane, by typical analysis. 

II. 3. Description of the Wind Tunnel 

The GAJLCIT hypersonic wind tunnel is a closed return, contin- 
uously operating tunnel with two test section legs. The legs are designed 
to be used alternately to provide for installation in one while the other 
is operating. Leg I lias nozzle blocks set for a nominal Mach number 
of 5. 8 in the 5" x 5. 25" x 29" test rhombus. The leg 1 operating 
limits are as follows: reservoir pressure, P Q = -5 to 100 psig; 

reservoir temperature, T Q = 225 to 325°F; for which Reynolds number 

4 4 

per inch. Re = 3. 88 x 10“ to 30. 3x10”. The Reynolds number per 
unit length is obtained by the formula 
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Re/L =<P 0 /> 0 ) M J { 



rf--rfz 



(T q /T) 



r-z 

?=T 



T 

H 
1 + 



198. 



to: 6 



.(la) 



T 



for a perfect gas, where zero subscripts refer to reservoir quantities, 
isentropic expansion from p is assumed, and the equation for the 
variation of viscosity v/ith temperature is the Sutherland formula: 



T ir 198. 6 

^2 ~ TTTW. 6" < T /V 



3/2 



For these tests with M = 5. 8 and T = 275 F: 

o 



( 2 ) 



Re/in = 2260 x P o (lb) 

where P q is measured in psia. The velocity of the test section flow at 
this nominal Mach number and is U = 2770 ft. /sec. as computed by 
the formula: 



<U/a/ 



2 

(2/M 2 ) CY- 1) 



» 




» 



where is the reservoir speed of sound . The reservoir pressure 
could be varied between approximately 0-100 psig, but critical 
starting and running conditions were avoided by limitations of 10 - 85 
psig. 



(3) 



A more complete description of the tunnel and compressor plant 
is to be found in References 19 and 20. A schematic drawing of the 
tunnel is included in Figure 1 and of the test section and related 
experimental equipment in Figure 2. 
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II. 4. Model Construction and Installation 

The basic model of the GALCIT hyper 3 oaic wake study at present 
is the two-dimensional circular cylinder. The body geometry was 
purposefully kept as elementary as possible in order to simplify the 
correlation of data. Once the complete flow field for this model is 
established the experimental methods for the simple n.odel can be 
extended and applied to more complex shapes. It was found that a 
0. 3" cylinder diameter was the maximum that could be employed for 
consistent starting of the tunnel flow. Thus, a . 3 M x 5 :l cylinder was 
mounted horizontally from wall to wall in the forward part of the 5 inch 
wide tunnel test section (Figure 2). Initially the model was installed at 
a point 13. Z" aft of the nozzle thi oat. In the later tests, however, the 
model was moved an additional 2. 6" aft to obtain better flow conditions. 
(See the discussion in Section III. ) 

McCarthy's total head and static pressure surveys were being 
made with a brass model. Although brass might be made porous, it 
was thought that a more suitable material should be found with a surface 
as smooth as possible and with natural random porosity. These proper- 
ties were found in an alumina refractory rod of the type used to support 
furnace heater elements. The refractory rods initially were of . 385 inch 
exterior diameter and 0. 13 inch interior diameter. The exterior of the 
model was carefully turned to a smooth, uniform . 3 inch diameter, but 
the interior diameter could be as much as . 003 inch non-concentric, 
because of warping of the rod. The slight non- symmetry of helium flow 
that would result from this interior eccentricity was considered of negli- 
gible concern as to its effect on wake measurements of gas concentration. 
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in view of other experimental difficulties that appeared to exclude such 
minor considerations at the time of model selection. Unfortunately, the 
rod material was somewhat brittle and the models had to be replaced 
frequently because of breakage. In spite of these deficiencies of the 
alumina rod the material was considered well suited to the present 
application, because it provided low, well distributed helium flow rates 
with low metering pressure during bench tests. Thus, the momentum of 
the ejected gas is small for a given mass flow, and should not appreci- 
ably alter the external flow and shock wave structure. It was desirable 
to achieve this result, if possible, in order to correlate diffusion data 
with the pressure data of McCarthy's experiments as directly as 
possible, and therefore obviate the necessity of duplicating a large 
amount cf his work. 

The bench tests of the models mentioned above included other 
considerations which are peculiar to this experiment. First, it was 
noted during preliminary tunnel runs that a light-surface coating of 
oil would collect on the model surface. The oil problem is a continuing 
difficulty in the GALCIT tunnel for which a satisfactory solution could 
not yet be found. It was therefore necessary to measure in some way 
the effect of the oil on the flow of helium through the rod surface. After 
several hours of running, the lightly oil-coated side of the rod was 
covered and the model was connected to a vacuum pump at one end, with 
the other end sealed. The difference in the amount of vacuum the pump 
would draw against a sealed space, and against the uncovered side of the 
model, was measured with a U-fcube manometer. The oil coated side 
showed a negligible change in porosity. The oil problem therefore was 
not considered serious, since the model could be rotated from test to 
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test to collect an even distribution of oil, and the flow level could be 
adjusted to provide the same flow of helium regardless of slight 
changes in model porosity. 

One must also make certain that all the rneiered helium is 
ejected through the cylinder and not lost in faulty seals at the tunnel 
walls. A clamping and locking arrangement and soft "o"-rings were 
installed to prevent such leakage. After each installation the vacuum 
test described above was applied to determine if the seal was con- 
sistently good. The helium was injected from both ends of the model in. 
order to provide a more uniform flow distribution across the cylinder 
span. The model was mounted between metal ports because of the 
difficulty of securing an injection model in glass ports when the in- 
stallation must be repeatedly disassembled. Schlieren studies were 
therefore not feasible with the cylinder injecting helium. However, 
comparison of the flow characteristics with and without helium flow 
could be obtained by taxing representative total head pressure traverses. 

XL 5. Control of Helium Injection 

The flow of helium entering the cylinder at the tunnel walls was 
metered with a standard Fisher- Porter Tri-flat flowrator, (Tube No. 
3-F-3/8-25-5/70, having a 3/S inch glass ball), in series combination 
with a smaller Fisher-Porter flowm.eter (Tube No. OIN-15, l/8 inch 
steel ball). The large flowrator provided stable measurement of the flow 
in spite of "back-pressure" variances of the model, while the small 
flowmeter provided a double check and a more sensitive, vernier -type 
monitoring of the flow quantity. The metering pressure was maintained 
by reference to a standard U-tube mercury manometer. (See Figure Z. ) 
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Although the temperature correction is small, the temperature of meter- 
ing was noted for each test. 

The metering of helium through the flowmeters was checked 
against the Fisher-Porter Company predictions and the previous 
laboratory calibrations for helium flow. This calibration was accom- 
plished by the timed displacement of a large volume of water at several 
helium flow rates. The large flowrator checked closely with the 
manufacturer's prediction, but the smaller flowrator was affected too 
seriously by back- pressure of the water to accomplish a good test. 

The flow of helium was therefore always based on the reading of the 
large flowmeter, with the smaller meter as a check against variances 
in cylinder porosity by the back pressure effect. The test was conclusive 
to about - 3. 5 per cent accuracy in the flow range of interest, which is 
also approximately the crder of accuracy in setting of the flow for the 
tunnel experiments. 

1L 6. Sampling Probe 

The probe employed to withdraw flow samples from the wake of 
the injection model was of the type designed for total pressure measure- 
ments in boundary layers. The tip was constructed by compressing a 
hypodermic tube to a flat orifice of 0. 004" depth with 0. 039" width. The 
tip was faired out to a 0. 25" C. D. stainless steel tube fitted with a 
0. 067" I. D. neoprene tubing at the downstream end for transport to the 
gas analysis cell. The same probe was used for total head traverses in 
the wake. 

The probe was mounted in a mechanism designed to traverse it 



axially from zero to 27 model diameters aft of the model centerline, and 
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vertically 1. 2" upward to 2. G M below the axial centerline of the tunnel. 

The probe could be set to an accuracy of about 0. Gl" axially and 0. 001” 
vertically. (See Figure 2. ) 

Samples were withdrawn despite the low pressure of the tunnel 
by the creation of a near vacuum condition in the conductivity cell. 

(This procedure will be explained in more detail in the next section. ) 

In the case of the lower tunnel pressures (near the model and at lov/er 
reservoir pressure) a relatively slow' time constant for gathering of the 
sample was noted - on the order of 30 seconds. It was considered 
essential to utilize the small probe, however, because of the narrow 
regions of helium distribution to be encountered near the model. 

II. 7. Thermal Conductivity Apparatus and Procedures 

The thermal conductivity or heat transfer capacity of a gas mixture 
is a distinct physical quantity that may be used in various ways to detect, 
identify and trace the concentration of one binary constituent in another. 

A very comprehensive summary of the principles and some of the more 

1 8 

common applications i3 presented by Baynes . An application in which 

j ^ 17 

helium was traced is the work of Rush, For stall, and Shapiro ' , 

v/here the method was utilized to trace diffusion in coaxial gas jets at 
low speed. The form of the equipment in the present experiment is 
similar to that designed by FvUsh and For stall. It was designed by 
M. D. Coffin at the C-ALCIT hypersonic laboratory. 

The thermal conductivity method of comparing two gas mixtures 
consists fundamentally in passing the mixtures in question through 
identical cell chambers with electrical resistance elements forming legs 
of a Wheatstone bridge. (See cell diagram. Figure 3, and wiring diagram. 
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Figure 4. ) The re ci stance elements in the two chambers have equal 

currents since the resistances are paired in the two legs of each 

chamber. Small variances can be removed by adjustment of the dividing 

resistance at the bridge input. When the gas conductivity in one chamber 

2 

equals that in the other, the wire heat, (I R), is conducted through the 

walls of the chamber, in equal amounts for the two chambers, and a 

balance is achieved with no voltage across the bridge. The resistance 

of the wires in the two chambers is not varied other than by the chamber 

temperatures. Also, in the bridge, one holds the current constant and 

essentially measures a change in E, the resistance of sample chamber 

wires. Thus R = RfT^ , where T w is the wire temperature. T is 

dependent on wire heat loss by the energy balance, 

• ... ~ 

V = <J g + % * I R . 

where neglecting radiation and in stagnant gas: 







heat loss through the gas = K ( at/6) 

T - , the difference between wire and chamber 

w c 

wail temperature 



5 = typical chamber dimension 

q Q = heat loss through wire ends = C -aT , where C = coefficient 
of wire conductivity times wire area, in appropriate units. 



Then 



K 



m 



i * * ( t- + c > < T w - T c> 



( 4 ) 



The chamber wails conduct rapidly enough to assume T = constant. 

c 



(In this case the chamber walls were constructed of brass. ) With con- 
stant current Eq. (4) then gives a relation between R, K m , and T w . 
The mixture conductivity, however, is a function of concentration of the 
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mixture species and the temperature, T , of the mixture, by the 

m 

kinetic theory relations (Cf. Reference 2, page 159). The mixture tem- 
perature of the combination of helium and air can be expressed as some 
function of a T, or as approximately the mean, T = {T w - T c )/2. One 
can thus replace K in Eq. (4) by the above equation from kinetic 
theory and then T w by a linear function of R, arriving at 

R = R q +aR a i{K Q + A K) . (5) 

Thus, each Wheatstone bridge resistance, R, is a function only of the 
change in concentration of the helium in an air-helium mixture surround- 
ing it. In this case, the bridge resistance in the cell cavity containing 
helium enriched air is reduced by the greater cooling capacity of the 
light gas. The variation of resistance in one chamber will cause a 
voltage unbalance across the bridge which can be read on a sensitive 
potentiometer. Equation (5) can best be evaluated in terms of this 
potential for purposes of calibration by laboratory measurement of 
known sample concentrations. {See Section II. 7. 1. ) Other schemes of 
operation are feasible, but this method is one of the most adaptable to 
this experiment and was especially convenient in that a commercially 
manufactured cell could be utilized. (The specifications and manu- 
facturer of the particular cell chosen are given in Figure 3. ) 

Although the fundamental procedures outlined above are con- 
ceptually simple, several additional problems arose in procuring 
representative samples from hypersonic flow fields, and in removing 
possible sources of error from the method. These are the considera- 
tions which determine specific design of the apparatus and operating 
procedure, and which will be discussed in succeeding paragraphs. 
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They can be classified under the following headings: 

(1) proper calibration of the cell 

(2) adequate handling of the sample 

(3) pressure differences betw'een the two chambers 

(4) temperature difference effects 

(5) instrument errors and accuracy. 

(An additional item of consideration here might be the handling of 
dissimilar mass flows between the two chambers. However, since 
hypersonic tunnel samples are widely variable in density and large 
errors would be introduced by convective cooling effects, it was decided 
from the outset to stagnate both sample and reference gases. ) The 
above five items will be discussed in turn with regard to how they affected 
the design and/or procedure of the experiment. 

II. 7. 1. Cell Calibration 

The calibration curve of potential across the bridge 
versus concentration of helium in a reference gas mu 3 t be obtained by a 
laboratory method, as shown in conjunction with £q. (5). The reference 
gas used must be one which remains consistent in heat conductivity after 
calibration for each subsequent filling of the reference chamber. In this 
case, if one were to employ room air in the reference chamber, filling 
the chamber for each period of use, an error may be introduced by 
departure from the calibration curve because of changes in the atmos- 
pheric mixture. The air mixture may change with small percentage 
variations of carbon dioxide and water vapor. The carbon dioxide 
variation will not be significant, but water vapor partial pressure is 
sufficiently variable and water vapor is sufficiently different in thermal 
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conductivity that noticeable error may be introduced if dry air is not 
used. (See Table I.) Since the air in the hypersonic tunnel is well 
dried and otherwise sufficiently uniform in content from one test period 
to the next, a sample of helium-free tunnel air was taken as the reference 
gas before each test. 

Having selected the reference gas as tunnel air, the cell 
calibration should be made with air-helium mixtures against dry air. 

The collection of large quantities of tunnel air is inconvenient, however, 
and the present calibration procedure involved the comparison of known 
mixtures of helium in dry nitrogen with pure nitrogen. Since the con- 
ductivity of nitrogen is very close to that of air, this procedure involves 
an error of smaller magnitude than that of preparation of the known mixtures 
and plotting and reading of the graph. The magnitude of the error involved 
can be estimated by noting that nitrogen is 0. 996 as conductive as air, 
(Table I). Thus, an error on the order of one part in 250, or . 4 per 
cent, is the result of comparing any sample against a reference of 
nitrogen instead of dry air. The difference in concentration obtained 
between comparison of a sample containing large quantities of nitrogen 
against a nitrogen reference and comparison of a similar air-helium 
sample against air is therefore quite negligible. 

The samples were prepared by allowing a bottle of high pressure 
nitrogen to discharge into a bottle of helium, with a sensitive gage 
measuring the pressure of the helium bottle accurately before and after 
the addition of the nitrogen. The partial pressures of the two gases in 
the volume then determined the concentration of the mixture. The 
preparation was therefore independent of possible nitrogen leakage 
during the transfer. The method of computing the concentrations is 
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based on the equation of state for components of a perfect gas mixture: 
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partial pressure of He, measured before mixing 

universal gas constant 
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He 



m 



He 



P He V m M He 

R n 

o m 



Similarly, 



m 



N- 



% V M N 2 
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Equation (6) is utilized to obtain the newly mixed concentration. 

A series of samples was prepared, having different concentrations of 
He in N-. As each sample was prepared, its concentration was read in 
terms of voltage deflection across the thermal conductivity bridge. The 
plot of these calibration points is a smooth curve with little scatter. 
(Figure 5). 



Temperature and pressure must be equal in the two cell cavities 
for a proper calibration. (Pleasure and temperature are discussed more 
completely in the sections to follow. ) Also the calibration curve is appli- 
cable only for the voltage and current utilized across the bridge when the 
sample deflections were read. The voltage and current choice depend on 
two conflicting considerations. Calibrations which are more sensitive 
and linear can be obtained for higher voltage-current combinations, as 
shown by the curves of Rush and Forstail 1 ^. Qn the other hand, it is 
desirable to keep cell resistance-wires warm by leaving the electric 
current on while drawing the next sample, in order to minimize the time 
required for the cell temperature to reach equilibrium in the measurements. 
However, during the sample collection, the gas density in the chamber 
becomes very low, and hence the voltage-current combination must be 
low so that the wires will not produce heat faster than it can be dissipated. 
(See Sections II. 6. and II. 7. 2. ) Even though the current was turned off 
before gathering a sample, the lowest sample densities are within a few 
millimeters of the vacuum pressure, when the current must be on to read 
the bridge deflection. It is therefore necessary to calibrate the cell with 
the lowest satisfactory power to prevent fusing of the hot wires. To 
accomplish the sensitivity adjustment in this case, a variable resistance 
across the potentiometer terminals is utilized. (See Figure 4. ) For 
calibration over the low concentration range of interest in the experiment, 
a current of 80 milli-amps with a 3 volt dry-cell circuit was employed, 
and the potentiometer resistance was set so that an essentially linear 
curve was obtained in the 0-2 per cent range. (See Figure 5 for equations 
of portions of the curve. ) 
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II. 7. 2. Sampling Equipment and Procedure 

The method of obtaining a representative gas sample from the 
hypersonic tunnel is somewhat more critical than in other applications 
because of the low pressures encountered (on the order of one millimeter 
of mercury). An intermittent sampling system is required in order to 
stagnate both sample and reference for proper comparison. 

As previously mentioned, the reference air was drawn from the 
tunnel in the present experiment prior to commencement of the helium 
flow. The bridge was balanced with the dry tunnel air in both sample 
and reference cavities. Samples were withdrawn, after the helium flow 
was commenced, by creating a near vacuum at the outlet of the cell 
(Valve 3, Figure 2). In order to compress the sample as much as 
possible to compare it properly with the reference air taken at atmos- 
pheric pressure, a hand-operated mercury pump was included between 
the tunnel and conductivity ceil. The pump in this case had a compress- 
ion ratio of about three to one. The pump included au-tube pressure 
scale, accurate to about 1 mm of mercury. Referring to Figure 2, 
valves i, 2, and 3 are sequentially opened and closed with a continuously 
operating vacuum pump so that the sampling steps are accomplished as 
follows: 

(1) evacuation of the cell, mercury pump, and connecting lines 

(2) washing of the new sample through the same volumes 

(3) re-evacuation of the volumes, which establishes a vacuum 
in the cell 

(4) drawing of sample with valve 3 closed, and the mercury 



flask lowered 
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(5) Compressing o£ the sample with valve 1 and 3 closed. The 
time required for each of the above steps is a function of the sample 
density. A low density sample enters the cell slowly because of low 
pressure differential between cell vacuum and probe total head. It is 
important to determine the proper time interval for each step. The time 
of taking one sample point is generally 3-4 minutes at beet by this 
procedure, but if it is hurried the sample may not be representative. 

The time for each step is found by taking an extended time at first, and 
then reducing the interval until a discrepancy is noticed in comparison 
with the same sample point taken over a longer time interval. This time 
interval is most critical for the lowest sample pressures. 

The sampling apparatus was thoroughly leak-checked prior to 
each test run by applying low vacuum to the mercury pump, tubing and 
cell cavities, closing the valves, and noting whether or not the vacuum 
is maintained by the mercury column. Even a slight leakage of room 
air into the system would cause gross contamination of the sample. Also, 
when a known sample is put in the sample cavity and left overnight, the 
same reading could be read again before the next test period, as a 
simple day to day check. 

Another pertinent remark should be made in regard to the method 
of collecting the reference sample. Since the pressure calibration 
depends on the difference between the cell cavities, one should be careful 
to see that the reference sample is collected at exactly atmospheric 
pressure. This objective may be accomplished by proper sequencing 
of the valves at the time of collection, (i. e. , shutting of the entrance 
valve an instant before the exit valve). Once the reference sample is 
collected at a known barometric pressure, it need not be changed again 
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during the test period. 



II. 7. 3. Pressure Difference Calibration 

The application of the thermal conductivity method of gas analysis 

to hypersonic flow requires some suitable means of correction for 

differences in pressure between the cample and reference cavities of 

the conductivity cell. For slight pressure differences this error is 

negligible, but when the reference gas is at atmospheric pressure and 

the sample at a few millimeters of mercury absolute, the error in 

neglecting the pressure difference may amount to 50 per cent or more 

of the concentration reading taken. Simple kinetic theory predicts that 

the heat conductivity coefficient of a gas is independent of pressure if 

the mean free path of the ga.3 molecules is much less than the chamber 

dimensions, and proportional to pressure in the opposite limiting case. 

The transition from one limiting case to the other is shown experimentally 
23 

by Bomelburg . His graph of heat loss versus Knudsen number 
(proportional to pressure) clearly indicates a smooth transition from the 
independent region to the proportional region at Knud sen numbers 
corresponding to the low pressures of the hypersonic wind tunnel, if 
the chamber dimensions are such as to permit such a transition. It 
was evident from the experiment by Bomelburg and from a few preliminary 
tests that the pressure correction to conductivity would be required in 
the present investigation. 

It is not simple, mechanically, to pump a representative sample 
from the wind tunnel at near vacuum to a pressure of one atmosphere in 
the sample chamber. Most methods of accomplishing this compression 
involve the use of a pumping fluid medium which may contaminate the 
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sample. While it is certainly desirable to investigate such means 

extensively, a simple expedient to circumvent the need for such a pump 

is to calibrate the error in bridge potential against the corresponding 

pressure difference. This calibration is accomplished by reading the 

bridge deflection of a known concentration sample at reduced pressures. 

The sample potential is initially checked at one atmosphere and then a 

portion removed in steps by a vacuum pump until very low pressures are 

reached. The new potential across the conductivity bridge is read for 

each step. A suitable means of plotting the resulting calibration curve 

for correction of tests is evident by physical reasoning, supported by 

23 

experiments like that of Bomelburg . One is led to the conclusion that 
the conductivity of the gas in each cell cavity, and therefore any unbalance 
of the bridge, is a function of the absolute pressure in the cavity, when 
the pressure is low enough that wail effects begin to become important. 
Thus, if 

p = pressure in sample cavity, cm. Hg. absolute 
s 

= barometric pressure of the open and of the U-tube column 
s measuring pressure in the sample cavity, cm. Hg. , 
a slowly varying function of time 



P 



a 



barometric pressure, p^ , of the room at the time when 
the reference gas is s collected, cm. Hg. 



4 p s height of mercury column, cm. 



Then 



P s = P b * 4 P 

S 

is the independent variable in a function which determines the emf error 

due to low pressure. Now p^ "= for any given calibration or test 

s r 

where the barometer does not change more than a few mm of Hg. , as 
is invariably the case. The emf correction is not detectable on the 
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instruments when the sample and reference cells are so nearly equal in 

pressure. One can thus say that p b = p._ so far as conductivity in the 

r s 

celt is concerned. To provide for barometer fluctuations, the absolute 
pressure of the sample cavity should be plotted in atmospheres as 

(P s /P b ) = (P b -Ap)/p b = 1 - (ap/p b ) 

s s s s 

for the abscissa of the calibration plot. The graph of Figure 6 is the 
result of the calibration. The points fail closely along the same curve 
for the three low concentrations utilised. Tunnel air reference points 
cannot be distinguished from nitrogen versus helium- nitrogen points in 
the small scatter at higher pressure differences. The scatter is the 
result of instrument inaccuracies, sampling variances, and the pressure 
measurement inaccuracy of the mercury pump scale. The pressure 
measurement was accurate to about 1 mm of mercury. Redesign and 
construction for better accux-acy was not considered essential in this 
initial wake survey, in view of other experimental inaccuracies of the 
same order of magnitude. Also, only a few points near the model 
involve pressure calibration at the very low pressures. A maximum of 
ten per cent inaccuracy for the pressure calibration at very low pressures 
was accepted in that it would normally involve much less than ten per cent 
error in the total readings of bridge potential. The corresponding per- 
centage error for any concentration level may be obtained from Figure 6. 
The scatter in the lower concentration region of the curve is entirely 
the result of slight variances in the sample and instruments used at the 
time of the calibration. The ernf correction here is small, and the best 
average of the points was taken in the curve fairing. (Later, the 
potentiometer was replaced and the micro-ammeter improved in order 



that less scatter would be obtained in taking the wake concentration 
surveys. See the instrument discussion in the next section. ) 

The pressure calibration certainly would 3 how a variation with 
concentration if the concentrations were of a higher magnitude than the 
3-5 per cent levels measured in this investigation. A series of pressure 
difference curves with various concentration levels as a parameter 
would be required for high concentrations. 

It is well to note that a pressure calibration for one geometry of 
the apparatus might not be the same as for another where the dimensions 
of the tubing and cavity were changed. This would be the result of 
mean free path and cell hot wire end loss effects. (These effects are 
discussed in the literature of kinetic theory and hot wire anemometry, 

C f. .Reference 22. } 

Finally, the following list of steps is given for a careful and 
precise pressure correction of the type that would be required for emf 
accuracies of 9" - 9S per cent with the lowest sample pressures to be 
found near the model at M = 5. S; 

(1) measurement of sample pressures in tiic range of 1 - 20 mm 
of m ercury absolute to an accuracy of 0. 1 mm, and higher 
pressures to 1 mm 

(2) accurate calibration curves for known sample concentrations 

(3) collection of reference gas at a known (atmospheric) pressure 

(4) monitoring of barometric pressure throughout the test 

(5) computation of pressure differences between cell cavities 
and normalization by the barometric pressure. (Since the 
variation in barometric pressure may be as much as 3 - 4 mm 
of mercury for a test run of several hours, each critical 
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reading can be computed using a curve of barometric pressure 
for the times of the test. ) 

The alternative to this procedure is the pumping of each sample 
to atmospheric pressure in the sample cavity with no contamination 
during the process. In the present test, the above five steps were applied 
as closely as possible, within the 1 mm accuracy of the mercury pump 
pressure scale, and with special care during the critical readings. 

IL 7. 4. Temperature Difference Effects 

Normal variances in* room temperature will not affect the 
accuracy of the thermal conductivity bridge as long as the temperature 
is equalized between the sample and reference cell cavities. (This was 
shown by Eq. (5) above, and its derivation. ) This statement can be made 
as long as the temperature coefficient of conduction for the gases to be 
compared is approximately equal. This factor determines the variation 
of heat conductivity with temperature, and is nearly the same for the 
gases considered here (Table I). The protection from radiating devices 
and stray air currents in the room by a wooden container and the use of 
high heat conductive brass for the cell material, help to insure the 
equalization of temperature. The brass also provides some equaliza- 
tion between cell cavities for samples brought in at a temperature 
unequal to the room temperature of the reference cell. The sample 
traverses a distance of 3 or 4 feet in the- room and is compressed in the 
mercury pump; also, there is a finite time between sample collection 
and reading of the bridge. The temperature of the average sample 
therefore cannot be greatly different from the room temperature of the 
reference cell. It is well, however, to eetimate the effect of a small 
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temperature difference in the bridge deflection reading. This estimate 

1 8 

has been accomplished by Daynes c , who states that the simple pre- 
cautions given above are quite adequate to prevent serious effect of 
temperature differences between ceil cavities. One may, however, 
install the cell in a thermo statically controlled housing if very accurate 
readings of low conductivity differentials are desired. 

Daynes, in his book, has given a very extensive coverage of the 
more obscure considerations in the method of thermal conductivity gas 
analysis. He includes such factors as vibration and other secondary 
effects. His conclusions regarding the detection of highly conductive 
gases in air, however, are that very small concentrations can be 
measured with accuracy if only the precautionary measures of the 
present experiment are employed. 

II. 7. 5. Control of Instrument Errors 

Referring to Figure 2 the instruments employed for measurement 
of the concentration potential difference were 

(1) the milli-ammeter 

(2) the thermal conductivity cell 

(3) the precision potentiometer with its associated standard 
cell and galvanometer (The galvanometer is utilized to 
obtain a null current through the potentiometer circuit. ). 

A simple method for control of most of the possible errors in these in- 
struments was available through the employment of the samples prepared 
for the calibration procedure. A sample v/as put in the sample cavity 
before each test run and the bridge emf was checked. The ratio of the 
calibrated value to the present emf was then a corrective factor for the 



readings of the bridge during the test. This procedure also served to 
correct for variations in room temperature, as a check on the bias 
setting of the bridge sensitivity, and as a check for damage or con- 
tamination of the conductivity cell. Along with a standard cell comparison, 
it also served to check the performance of the potentiometer. Carefully 
mixed samples of known concentration are therefore indispensable in the 
use of the thermal conductivity method. Frequently, as mentioned in the 
sample leakage control discussion above, the known sample was left in 
the cell for extended periods of time and the same reading was obtained 
upon recheck. The reproducibility of such readings on known samples 
was within about one per cent at all times. 

The bridge was balanced by a null indication on the galvanometer 
and potentiometer before each test with the dividing resistance at the 
bridge input, (B+, Figure 4). This balance was accomplished with 
nitrogen in reference and sample cavities, and then with air in both as a 
recheck while gathering the tunnel reference air before injection of 
helium. The stability of the ceil zero potential was within 0. 05 m. v. 
at all times. 

XI. 8. Total Pressure Apparatus and Procedures 

The magnitudes of injected helium concentration to be measured 
in a flow depend inversely on the local density and velocity of the flow. 

The determination of these quantities at any station of interest is 
therefore mandatory to a complete understanding of the diffusion 
processes. It is not an easy matter to obtain such quantities in a 
compressible flow . The total and static pressures are relatively 
easy to obtain, but the measurement of another state variable is 



required to determine density and velocity. The additional measurements 

might be the speed of sound, the temperature, or the direct measurement 

of velocity and density, if these measurements were feasible and 

dependable. At present, methods for close determination of these 

flow variables are not fully reliable. One may assume isoenergetic 

flow, however, without serious error in first order calculations, and 

compute an approximate local temperature. The pressures are then 

measured with suitably designed probes and associated manometers 

or transducers, after which the Mach number can be obtained through 
ZZ 

the Rayleigh pitot relation 



Y = ratio of specific heats. 

The Mach number and pressures may then be employed to determine 
other fluid quantities. 

As previously mentioned, a pressure study by J. McCarthy in 
the wake of a 0. 3 inch brass cylinder was also in progress at the time 
of the present experiment. It was desirable in terms of wind tunnel 
running time to utilize the pressure data from the brass model in the 
diffusion studies made with the porous model. Inis procedure required 
a comparison of the two models as to flow characteristics, without 




$ 



( 7 ) 



where 



p = static pressure 



p^ = total pressure behind probe normal shock wave 

z 

= Mach number ahead of probe shock wave 



29 



injection and with helium injection at the level to be metered during 
the wake diffusion tests. The total pressure profiles were considered 
the best means of comparison. Total pressure profiles behind the 
injection model were therefore made with the equipment employed by 
McCarthy. This equipment will be described more completely in his 
report, but a brief description will be given here. 

The probe utilized in the total head traverses behind the cylinders 
was the same instrument employed to gather diffusion samples, and was 
described in Section II. 6. The total pressure in the probe was first read 
with a mercury micro-manometer at one or two points in the flow, after 
the tunnel had been stabilized at the test reservoir conditions. The total 
pressure was then diverted to a Statham 1 psi differential pressure trans- 
ducer with a reference pressure provided by a silicon oil column. The 
transducer converted the pressure to electrical potential which was 
traced on a Mosely xy- plotter. The transducer-plotter system was 
calibrated by measuring the pressure with the mercury micromancmeter 
at a few points. The output from a voltage divider coupled to the probe 
traversing mechanism was fed into the x- scale of the plotter. Sufficient 
additional micro-manometer points were read to establish that the 
recorder scale had not shifted during the tunnel traverse. All pressures 
were referred to a very low vacuum and therefore were plotted as absolute 
pressures. The micromanometer pressures were non-dimensionalized 
with the tunnel reservoir pressure converted to absolute units before 
making the traverse. Thus, the pressure traverses are plotted in a 
convenient form for comparison of flow characteristics. 
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IvlcCarthy had also taken static pressure traverses of the wake 
behind the brass cylinder. A comparison of static pressures with the 
porous model, however, was considered unnecessary in tliat total head 
profiles would give sufficient insight into any flow differences, which 
should not be large. 

Results of the wake studies depend greatly on whether two 
dimensional flow actually existed in the tunnel across the cylinder span. 

A few representative total head traverses were therefore made to check 
the two dimensionality of the flow approximately one inch from the tunnel 
centerline. 

II. 9» Correction and Correlation of Bata 

The reduction of emf data from the conductivity bridge requires 
a certain sequence of steps according to the discussion concerning the 
several items of equipment in the preceding sections. Before outlining 
the reduction technique, however, it will be necessary to show how the 
build-up of helium: in the tunnel circuit may affect the data. 

The quantity of helium in the tunnel air can be predicted by a 
brief analysis. Helium concentration does not build up in value indef- 
initely, but approaches an asymptotic ratio of helium mass to tunnel air 
mass which is dependent on the amount of make-up air added to the tunnel 
circuit and the corresponding leakage rate from the circuit. The 
governing equations are analogous to those for water vapor in the tunnel 
air, where the amount of water vapor is controlled by the bypassing of 
circuit air for drying to prevent condensation in the test section. The 
equation for the tunnel helium concentration can be derived as follows: 
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Let m^ = constant mass/time of He entering the tunnel 

= constant mass/ time of make-up air entering the tunnel 
= mass of mixture in the tunnel, 

= concentration by mass of He in the mixture 



M 

K 

rli 



= mass of mixture leaving the tunnel with concentration R, 
which must equal the total mass entering for steady tunnel 
operation, or 

rfj = m_ T + m = a constant 
He a 

Therefore, the rate of change of helium concentration is given by 



m. 



... - {m + m ) K 

dx/dt = —£52 He a 



sr 



The solution of this equation with R = 0 at t = 0 is as follows: 
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He 



ni He ■ tt a 
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1 -exp{- ne a t) 

M 
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Equation (8) demonstrates that the tunnel concentration tends to an 
asymptotic constant value as t becomes sufficiently large. This result 
was confirmed experimentally in preliminary teats. Within the accuracy 
of the thermal conductivity measurement, the time required to reach a 
constant value was on the order of 10 - 15 minutes for the normal quantity 
of make-up air, rh , utilized to replenish the tunnel circuit. The tunnel 
concentration with normal make-up air and low reservoir pressure is of 
the same order of magnitude as the lower concentration readings. It is 
therefore desirable that as much as possible of the background concen- 
tration be removed. Make-up air was increased to accomplish this 
removal by excess leakage until levels of helium concentration about 
. 06 per cent or less were obtained at the several reservoir pressures. 

Equation { 8 ) shows that the asymptotic constant value of tunnel 
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concentration reached is independent of tunnel P q , and is dependent only 
on the rate of helium and air added to the tunnel. (The rapidity in reach* 
ing the asymptotic value, however, is dependent on through the tunnel 
cycle mass of air, M. ) This result was also shown to be correct when 
the valves admitting make-up air were left fixed and the tunnel was 
changed for some of the tests. IV hen the mas8 rate of the helium and the 
corresponding tunnel concentration are known, the make-up air rate can 
be computed by Eq. (3). If the make-up air rate is too high the tunnel 
air will not be sufficiently dried, and errors in measurement of helium 
concentration may be introduced through lowering of the air conductivity 
by water vapor. The tunnel concentration may vary during a test if the 
make-up air is changed. It was therefore checked at intervals and 
corrections were made as necessary in the data reduction. 

Perhaps the best method of discussing the reduction of concentra- 
tion test data is to follow a typical data sheet (Table II) and reduction sheet 

(Table III). The values of concentration in millivolts read for a series of 

t 

distances, y, from the vertical centerline are entered in the reduction 
sheet, along with the corresponding pressure differences, a p, between the 
current atmospheric pressure and the sample cavity pressure. The com- 
putation is started in column 6 with the subtraction of the tunnel background 
concentration from the emf reading. The tunnel concentration sample is 
always compressed to atmospheric pressure, which allows subtraction 
before a p correction. The computation is continued with column 7, where 
the emf reading of the conductivity bridge is corrected for instrument 
errors by the ratio of the emf of a known sample to its current reading. 

(The instrument correction is not large, for it can be controlled by the 
variable resistance across the potentiometer terminals. ) This ratio 
corrects for variations in current measurement, room temperature. 
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and slight change in cell characteristics as discussed in previous 
sections. 

Column 9 of Table III scales the a p vacuum versus barometric 

pressure reading of the mercury column to that barometric pressure. 

The barometric pressure taken al one time during the experiment is 

sufficiently accurate unless very low tunnel densities are expected. If 

these densities are such that the ratio ip/p. > . 99 , it is well to take 

bs 

into account shifts in barometric pressure of greater than 1 mm. 

Fortunately such low pressures are found only very close to the model 
in restricted regions, where the concentration is normally high. The 
percentage error is therefore low in the total emf reading for a relatively 
high error in a emf caused by cell cavity pressure differences. 

The a, p/p, „ value is utilized in Figure 6 to obtain the £>enif 
correction for column 10 of Table III. Finally the emf versus concentration 
graph of Figure 5 is entered to obtain the concentration, I', in per cent 
of helium by mass. 

The remaining step before plotting of the concentration data is 
the determination of the centerline from the maximum concentration 
point. This step is carried out after the a p correction lias been applied 
in order to determine the true maximum. In some of the tests only 
enough points were taken vertically above and below the centerline to 
determine the true maximum for rate of decay versus axial downstream 
distance. If a profile were somewhat skewed, the centerline was 
determined by equal spacing of half-maximum points. (This point will 
be covered in more detail in the later discussion. ) 

The helium flow for some of the early tests was not exactly the 
0. 0031 Ib/min of the later tests. It was changed after the profile runs 
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to improve the exactness of the setting. A scale factor of 1. 15 was 
applied to the data of the early runs to refer all runs to the same helium 
flow. The linearity of such a correction for small ranges is shown by 
Figure 7, which is a plot of a few points taken at different helium flow 
rates with the probe fixed at an arbitrary location behind the model. 

This procedure would be in error if the concentrations measured were 
not in the linear portion of the millivolt versus per cent calibration curve. 
Figure 5. 

Finally, before proceeding to the results of the experiment, an 
estimate should be made of the errors to be expected in the concentration 
data produced by the equipment and procedures employed. The types 
and causes of the errors have been covered to a great extent in the above 
discussion. A more complete summary of them is included in Table IV, 
with an estimate of magnitudes wherever possible. This table may serve 
as a guide to improvement of the accuracy of the method. A fairly com- 
plete reference concerning accuracy of thermal conductivity measurement 
is the book by Daynes^. 

IL 10. Summary of Test Parameter Variations and Test Objectives 

In view of the simplicity of the apparatus in the present exploratory 
investigation and the resultant long time intervals required in the 
measurement of individual wake samples, the number of parameters to 
be varied in the tests was held to a minimum consistent with the basic 
objectives previously cited. The most important parameter in the wake 
flow processes of interest is of course the Reynolds number. Variation 
of Reynolds number was achieved by holding the reservoir temperature 
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at T q s 275°F and altering the tunnel reservoir pressure alone, 
according to Eq. lb. (Reynolds numbers of the tests are included in 
Table V. ) The helium ejection rate was held at a constant low nominal 
value of . 0031 lb/min throughout the tests. This was ascertained to be 
the flow quantity desirable for proper tracing of the processes with 
minimum effect on the momentum of the basic flow. 

The investigation, with the above considerations, included the 
following tests: 

(1) Determination of the maximum concentrations behind the 
model for representative axial distances to the limits of probe travel 
(27 diameters) for supply pressures P Q = 10, 35, 60, 85 psig, and all 
other conditions held constant. The maximum concentration would be 
determined by taking three or four samples near the tunnel centerline 
until the maximum is well defined. These tests were required to 
determine how the maximum concentration decays with axial distance 
behind the model. 

(2) Measurement of flow samples for concentration at a sufficient 
number of points off the vertical centerline to determine a concentration 
profile at various downstream stations. These traverses for a number of 
such stations would then display conditions of spreading and regions of 
similarity in the flow. The combination of centerline maximum con- 
centration points taken and the profiles should then determine regions of 
turbulent flow. To restrict the amount of tunnel operating time, the profile 
traverses were limited to the extreme pressures considered, P Q = 10 and 
85 psig. 

(3) Total pressure traverses were planned at a representative 
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downstream distance, x/D * 5, for the supply pressures F q = 10, 35, 60, 
and 85 psig, with and without helium injection. Comparison with data 
from the brass model would then show corrections that should be made 
for pressure differences, in utilization of that data for the porous model. 

(4) Total pressure traverses were to be taken off centerline to 
determine if the flow is two- dimensional at a representative station, 
x/D = 5 for P q s 10 and 85 psig. 

(5) Monitoring of the flow at intervals during the concentration 
runs was also to be included, by the taking of sufficient total head 
traverses to see if conditions in the flow field remained consistent. 

These were considered the minimal tests that could be made to 
determine a satisfactory general description of the flow and diffusion 
processes affecting the injected helium. From such a limited investi- 
gation, then, a much more extensive survey could be planned, with 
improved equipment and technique for coverage in reduced tunnel 
operating time. 

Each te 3 t was designed to duplicate points of a previous test to 
check for repeatability of helium flow setting, model installation character- 
istics, conductivity apparatus variations, and possible alteration of the 
basic flow. Correlation of such repeated data would then establish a 
confidence level and accuracy check of the procedures. 

A typical data sheet for the recording of concentration test data 
is included as Table II. The fixed parameters are included in the sheet, 
and the remaining blanks are explained by the discussion of the previous 



sections. 
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III. RESULTS AND DISCUSSION 
III. 1. General Discussion of tho I ’low 

At hypersonic speeds the flow field around a blunt body is domin- 
ated by the bow shock (Figures 8, 9). The static pressure on the front 
portion of a circular cylinder with its axis perpendicular to such a stream 
is much greater than ambient pressure, and the pressure falls steadily 
along the surface. When the boundary layer is laminar, flow separation 
occurs son e dietance aft of the position of maximum thickness (Figures 
Q, 9). The static pressure along the inclined free shear layer is very 
nearly constant and this layer is nearly straight until just upstream of 
the characteristic "neck" formed in the wake. Then the flow deflects 
and the pressure rises abruptly in a short distance, while the flow velocity 
in the wake cavity along the dividing "zero" streamline is brought to rest. 

A second shock wave is produced by this flow deflection. This oblique 
wave is intercepted by the expansion fan generated at the body surface 
and decays rapidly in the downstream direction. 

When the free shear layer is laminar the width of the "neck" is of 
the same order as the boundary layer thickness at the separation point 
on the body. In the Reynolds number range of the present experiment, 
the cylinder boundary layor io laminar at separation, which is confirmed 
by the fact that the neck is exltcn cly narrow in con pariaon to the body 
diameter. Rapid lateral diffusion of vorticity and heat takes place as ttti 
fluid moves downstream. Of course, the gradients of enthalpy and velocity 
in the "external" flow which has traversed the bow shock are also being 
smoothed out by laminar diffusion, but the time scale of this process is 
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several orders of magnitude longer. 

If a foreign gas is ejected from the surface of the circular cylinder 
it will also be concentrated at the narrow "neck", provided the free shear 
layer is laminar. Downstream of the neck, the lateral diffusion- of this 
tracer gas furnishes a measure of the efficiency of the laminar mixing 
process. At sufficiently high Reynolds numbers, transition to turbulent 
flow is observed to occur in this inner wake downstream of the neck. 

The subsequent lateral diffusion of the tracer gas in this case gives 3 ome 
insight into the rate of spreading of the turbulence into the "external" 
flow field. Of course the problem is complicated near the neck, but 
further downstream the lateral concentration profiles should exhibit a 
certain similarity. Consider a point far enough downstream of the neck so 
that u - u < < u . Then the Oseen approximation can be employed, and 
the diffusion equation to this approximation is as follows: 



where D T , = , the binary diffusion coefficient for laminar flow, and 

D r , = D,-, , the turbulent eddy diffusion coefficient for turbulent flow. For 
laminar flowaie absorbs the density into the Howarth-Dorodnitzyn variable 
by lettins 



pu e (8K/8x) = (8/eyH p ) 



(9) 




o 



By utilizing the approximation 



2 

pip D l^) = constant = C p £ (£>ig) e 



as analogous to the usual py< = constant for the velocity defect solution, 
Eq. (9) is reduced to the form: 
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3K/8x = 



C ( D 12>e 



8^'K 



u 



(10) 



e BY 

In addition there is the condition for the conservation of the tracer gas, 
namely, 

r° 

n, = \ puK dy = constant, (11a) 



or 
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K dY = constant 



cub) 



The solution to Eq. (10) with the condition given by Eqs. (11a) and (lib) 
is the well-known Gaussian distribution 

.2 



K = 



/cr 



e * p 1 ' ’ • (12) 



where 



i. e. , 



A 



= Aj (m/p e u e ) , (Aj = a constant) ; 



X <= A, 



m 



exp ( 



C p (D, -j) u x 
' 12'e e 



4(D u ) e C(Vu e ) 

The conclusions indicated by this first order laminar solution are that 
< a > K n : a,:“ l/ ^ , given p e , rh ; 

(b) Since (D u ) e ~ ^e^/e/Pe) . 

K max^ l/ ft e ’ ^ iven ™ 1 



) .(13) 



(<=) K /^ m ax - 8X ? l “ ^[L» 12 / e <3(x/u e ) 



) 



indicating the Gaussian shape of the concentration profile. 
Now for turbulent flow the situation is more complicated because 
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the turbulent eddy diffusion coefficient is not independent of axial distance. 
Roughly speaking, D™ = Bau 5 rr , whereau is the velocity difference 
across the turbulent inner wake, B = a constant, and S T , is the breadth 
of the inner wake. As the turbulent fluid penetrates into the "external flow", 
Drj- increases. Nevertheless a first approximation for the turbulent case 
is obtained by regarding D,-, as a constant. Making another first order 
approximation, let p = , a mean value in Eq. (9). Then the lateral 

concentration profile for the turbulent case is similar to the laminar 
solution, except that ) g is replaced by , Y is replaced by y, 
and C - 1 in Eq. (13). Nov/ 




B(C-) D 

' D w 

“3 



> 



v/here (C-„) i9 the drag coefficient associated with the momentum defect 
of the turbulent inner wake and D is the cylinder diameter. Then for 
turbulent flow; 



K = A- 



m 

“p u IT 
*e e 



I 



B C 
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— exp ( 

/x7S B C D 



Dx 



) . (14) 
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The conclusions drawn from the first order "solution" of the turbulent 
case are that 



( a ) \ ax ^ 1/ /x > given p c , rh , ; 

w 

< b > K max~ */Pe > given x/D , m , 

w 

(c) The profile of K/K , is Gaussian. 



Before examining the experimental data in the light of the 
preceding rough similarity considerations, certain flow conditions 
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discovered during the tests must first be discussed. 

IIL 2. Tunnel Flow Conditions and Variations 

During the course of the pressure traverse tests made with the 
brass model and the concentration tests with the porous model, certain 
discrepancies began to appear. (Both models were mounted in the same 
tunnel position, 13. 2 inches from the nozzle throat. ) Discrepancies 
were first seen in the pressure data taken by McCarthy, prior to the 
concentration runs. (See Table VL ) In particular, it was noted that 
reproducibility of the pressure traces was not good, and that square-shaped 
and non- symmetrical wakes were being obtained. Actually, Gaussian 
distributions of wake diffusion quantities are to be expected, according 
to the discussion in Section III. 1. The discrepancies in the pressure 
traverses were attributed in part to oil droplet formation on the brass 
model. (It diffused over the surface of the porous model. ) However, 
it was not until a quantity of pressure and concentration data was taken 
that the discrepancy affected the data seriously enough to be certain that 
it, in fact, existed. The chronological history of the flow can be followed 
by referring to Tables VI and VII and the corresponding figures (Figures 
10-21), which are the complete presentation of the experimental data 
obtained. In particular. Figure 13 shows the effect of the flow change 
in spreading of helium mass; Figure 14 shows the effect on axial decay; 
and Figures 16-18 are pressure traverses which show the variations in 
the momentum defect of the wake core. Figure 19 demonstrates that the 
flow is identical for the two models. Figures 18 and 19 demonstrate that 
the effect of helium flow is that of a slight lowering of the wake total 



pressure, by 5 per cent or less at P q = 85, and widening of the shocks 
and wake core at P q =10 psig. Density effects of the light gas are 
undoubtedly the factors which caused these differences in total head 
traverses. 

When all of the above figures are reviewed, one comes to the 
conclusion that tunnel flow perturbation is certainly evident. A few 
pressure traverses were taken without the model. These traverses 
and some calibration history of the tunnel revealed that two small waves 
originating in the nozzle throat region crossed the tunnel at almost the 
exact position of the model. Their symmetry had made them impossible 
to detect in the early data. Their effect was to cause intermittent 
variation in the pressure of the wake and widening of all the wake total 
pressure "buckets". This effect can be verified by noting that the 
pressure traverses taken with the probe off center actually correspond 
closely to those taken with the model moved rearward Z. 6 in. out of the 
perturbation region. (See Figures 16c, 18b, 18c, and 19b. ) These con- 
clusions were confirmed by intermittently taking concentration data and 
pressure traverses with the model at the rear position until the data of 
concentration runs 10-15 was obtained without another shift of the flow. 
The runs taken with a stable flow then serve as a means of correlating 
the previous concentration profile data and obtaining some useful 
information from it. This procedure is justified only because of the 
precautions taken in all concentration tests, namely those of repeating 
sample points from preceding runs to check reproducibility. The small 
variations in the duplicated points which were noted at the time of taking 
the profile data could now be understood. Originally, these were attri- 
buted to errors in helium flow setting or to instrument errors. 
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Ii was considered impractical in this exploratory study to repeat 
the data of the first runs for the present investigation, after moving the 
model to the stable position. Rather, it was decided to utilize the data 
in its present form to extract as much information as possible, and then 
to repeat and expand the experiment in another investigation, with im- 
proved equipment based on the knowledge gained to date. Actually the 
effect of this compromise on the remaining discussion is not as great 
as might at first be expected. 



III. 3. Continuity Check of Helium Flow 



By integrating the diffusion equation to suitable limits of the 
concentration profile, one obtains 





+ 6 

(8/8x)(paK)dy + 

*°K 



+o K 

(8/ 8y)(puK)dy = f (8/8y)(pD 12 Tj~)dy • 

“ 6 k 



Since the concentration and concentration gradient in the y- direction are 
both zero at the wake edges, the integral equation is reduced to 




puK dy 



= 0 



> 



where the order of differentiation and integration has been interchanged. 
This result can only be valid if 




constant = rh^/ unit length 



(15) 



which expresses the fact that the mass of helium injected from the cylinder 
into the stream must equal the integral of the puK profile taken at any 
station, x , and integrated over the span. Thus if one assumes two- 
dimensionality of the flow, the integral 
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puK dy dz = rh^ , mass/sec (16) 

° - 6 k 

would e^ual the helium flow quantity metered into the cylinder. The 
integral provides a means of checking whether or not the concentration 
measurement is accurate, providing the product pu can be determined. 

The method of computation for pu follows from the equation of state and 
appropriate modifications: 




pu = p (Ma/RT) = pM ; (17) 

where Y m and R m are the ratio of specific heats and gas constant of the 
mixture of helium and air. If the assumption is made that the flow is 
iso-energetic (with possible error not greater than approximately 10 per 
cent), the local temperature can be replaced by the well-known relation. 



T = 



1 + 



T 

o 

t~=t 

m 



(18) 



M 



The Mach number may be determined from the pressure data taken by 
McCarthy, which corresponds to the concentration profiles taken before 
large shifts in the tunnel flow. The method of obtaining the Mach number 
from this data is given in Section II. 8, for Y - V . Although the static 
pressure is given only for the wake center, McCarthy has found by taking 
representative traverses that static pressure is very nearly constant over 
wake cross-sections between inner shocks downstream of the wake "neck". 

To provide greater accuracy in the puK calculation, one should 



apply a correction to the static pressure for the helium defect in the wake. 
(See Figures 18 and 19.) A measurement of the static pressure with and 
without helium flow should ordinarily have been made. It is reasonable. 
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however, that the velocity is not changed and that the 3tatic pressure, p, 

g 

is changed by the same ratio as the total pressure, P, = p + p/ U , 

i? 2 

because of the density effect of the helium. The effect in this case was 
so slight that no correction to p was made. The Mach number does not 
change in any case, with the assumption that the velocity change is very 
small. 



An accurate computation of the puK quantity also requires that 

the parameters ^m and R be replaced by their true values, rather than 

m 

by approximating them with Va = 1.4 and R = 1, 716 ft. "/ sec 4 °R for air, 

cL 

because the parameters of the two binary constituents are quite different in 
the case of helium and air. (R u ._ =12, 438 ft. ^/sec^ °R and = 1. 667) 

The equation for the gas constant is 
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where K = concentration of helium in air, by mass. The equation for the 
ratio of specific heats is c 
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1 - K + K 

When the numerical values for the gas parameters are substituted, the 
two equations become 



R m = 


1, 716 + 107. 22K 




(19b) 


II 


140 + 5. 351K 
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• 


(20b) 



The value of puK obtained for each vertical distance, y, by the 
above procedure, is plotted as in Figure 23. The integration by planimeter 
then gives the mass per unit span. The double integration is complete when 
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the mass/ spaa is multiplied by the span length, b = 5 in. The two 
example curves of Figure 22 came from the profile data of Figures 10 
and 11 and the corresponding pressure and Mach number profiles of 
Figures 20 and 21. The non-symmetry of the puX profiles is to a certain 
extent the result of non- symmetry in the total pressure profiles, but 
may also be caused by the non-uniform injection of helium through the 
cylinder. Four such integrations were possible with the available data 
and were completed with the result: 



P o = 


85 


x/D = 


15 


rfuj e = . 00279 lb/min 


o’ 0 

II 


85 


x/D = 


20 


m_, = . 00300 lb/min 

n6 


P o = 


10 


x/D = 


9 


HiRg - • 00307 lb/min 


P o = 


10 


< 

t) 

II 


15 


m.^ e s . 00281 lb/min 



These values are to be compared with the flowrator setting m^ = • 00305 
- . 00010 lb/min. (The second value above was taken by doubling a half- 
concentration profile, after the flow had shifted slightly. ) Even though 
the several approximations above were made, the values have less than 
10 per cent error and 10 per cent spread. It could be argued that the 
results are somewhat fortuitous, but it is unlikely that gross errors 
could have cancelled each other for the four different cases at two tunnel 
reservoir pressures. One is led to the fairly safe conclusion that: 

(a) The iso-energetic temperature assumed is within about 10 
per cent of the true value. 

(b) The static and total pressure measurements are dependable. 

(c) The concentration measurements are also to be trusted 
within the accuracy predicted by Table IV. 

(d) The tunnel instability was not serious in this case because 



47 



the pressure and concentration profiles corresponded. 

(e) Even better accuracy should be obtained if the instrumentation 
is improved and the model is located in the rear position. 

(0 Two dimensional flow is, for all practical purposes, a 
certainty. (This conclusion is supported by the relatively good agreement 
between the forward off-center pressure traverses and those taken with 
the model in the rear position. ) 

III. 4. Discussion of the Diffusion and Similarity Properties 

The values of all maximum concentrations taken at the wake 

centerline are plotted in Figures 14 and 15. The plot of Figure 14 was 

included to show the effects of the flow shift on axial decay at P =85 psig, 

o 

when the model was in the forward position. The points plotted in Figure 
15, however, are those taken with stable flow conditions when the model 
was moved to the rear position. Unfortunately this shift in model position 
reduced the probe travel behind the model centerline to 27 diameters 
instead of the previous 36, but by this expedient a complete and con- 
sistent set of axial decay points at four P values was obtained. The 
flow was checked before, during, and after the runs by total pressure 
traverses with the probe at 5 diameters. The stable traverses of Figures 
16 - 19 were consistently duplicated. Having established a fairly high 
degree of confidence in the data of Figure 15 by the continuity check of 
the previous section and the evidence of stable flow, this plot may be 
discussed in the light of previous experiment and theory. 

When the word similarity is used in the following discussion it 
should be construed to mean similarity or approximate similarity in the 
general sense; that is, a scaling of the flow according to certain laws 
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relating the coordinates and diffusion parameters, such that all cross- 
section profiles taken at downstream positions have the same structure 
and shape and follow the same decay rule. Townsend attached this 
meaning to the term "self-preserving"^ , and reserved the term 
"similarity" for the more specific Reynolds number similarity in 
turbulent flow, requiring that the "processes which determine the main 
structure of turbulent motion are substantially independent of viscosity". 
This essentially applies to all turbulent free boundary flows in the con- 
tinuum regimes at high velocity. It does not mean, however, that the 
boundary conditions of the flow cannot be influenced by viscosity, as 
indeed they must in the case of the cylinder wake, where the drag 
coefficient is a function of Reynolds number. 

III. 4. 1, Downstream Decay of Maximum Concentration 

The theory shows for the present case that the maximum helium 
concentration should be distributed down the tunnel centerline in inverse 
proportion to the square root of x , the axial coordinate. Further 
K „ ^ x 3 should be true of the flow in either the laminar or turbulent 
case. The plot of Figure 15 indicates that this condition is true in only 
the P q s 10 psig case for the region beyond x/D =4. A constant slope 
for the other reservoir pressures is indicated, however, for regions 
sufficiently far downstream. In the physical case, it is realized that the 
model axis is not necessarily the origin of the similar flow, and that in * 
order to compare the slope of the decay curve with theory, one must 
determine the virtual origin of the similar flow by a suitable auxiliary 
plot. This procedure is accomplished for the present data by plotting 
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the values l/K ^ against the coordinate, x , on a rectangular plot. 

(See example. Figure 23. ) Though scatter is enlarged by this procedure, 
a straight line through the points is evidently a good approximation to 
their distribution, showing that the x" 2 condition is satisfactorily met, 
with virtual origin as indicated by the x-intercept. The resulting 
correction to the plot of Figure 15 for all four reservoir pressures is 
included in Figure 24. It can be seen that the maximum concentration 
does at least closely approximate the decay condition given by the theory. 

It should be noted that a few points in Figure 15 fall below the 
normal scatter of the data. A close check of the experimental data sheets 
revealed that the points m.arked in the figure were taken at the end of a 
day of operation with high tunnel make-up air. It is quite likely that the 
tunnel drying equipment could not remove sufficient quantities of water 
vapor in these cases, and the conductivity of the tunnel air was lowered 
as compared to the reference sample taken at the beginning of the run. 

With the relative conductivity of water vapor as given in Table I, about 
1 per cent by mass would have been sufficient to cause the observed 
deviations. This is then a likely reason for the appearance of the profile 
taken at 24 diameters in the =10 psig case of the early data (Figure 11) 
which was also taken at the end of a long test period with high make-up 
air. The use of high make-up air for other than short runs therefore 
appears questionable. Another expedient would be to collect a reference 
sample for each tunnel sample, but this would increase the time of 
taking a reading and the possibility of error in the procedure. 
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ILL 4, 2. Variation of Maximum Concentration with 
Tunnel Pressure 

The separation of all the curves in Figure 24, for the various 
reservoir pressures, does not appear to be a function of alone, but 
is dependent on some other condition of the flow. Since the theory calls 
for 1/ /p” proportionality at a given x/D in the laminar case, one 
immediately checks to see if the condition is satisfied in the present data. 
This is done by checking the slope of the auxilliary plot made to determine 
the virtual origin, x q , or by checking the ratio between straight-line 
portions of the curves in Figure 24. The condition is satisfied for the 
case of P Q = 10 versus P Q = 35 psig, with 

rP o ('35)/r o (10) = y(35 + 14.4)7(10+14.41 = 1.42 , 

where 14. 4 psi is the average barometric pressure of the laboratory 
elevation. It is not satisfied for 

or . 

where the values for the absolute square root ratio are 1. 23 and 1.16, 
respectively. Furthermore, the flow does not satisfy the predicted 
condition for the fully turbulent case in either of these comparisons, 
in that the ax values arc not scaled directly by the pressure ratios, 

P o <60)/ p o <35) and P o (85)/ p o (60) 

A possible conclusion indicated by these facts is that the flow at 
P q s 60 psig is in a region of transition from laminar to turbulent flow 
and that the flow at P q = 85 psig is fully turbulent. The flows at P^ = 10 
and 35 psig are of sufficiently low Reynolds number that they are laminar. 



These statements arc supported by cue related pressure traverses which 
were taken intermittently with, the concentration runs to check flow 
stability (Figures 16 - 18, runs 10 - 13). The wake portion of the pressure 
traverse in the P q = 60 psig case is narrower and deeper than that of the 

P = 85 trace. This trend is opposite to that between the P =10 and 

o o 

P q = 35 psig cases, where the P Q =10 trace becomes shallower at the 

wake center on the nondimensional scale, and remains at abou: Lhe same 

width. Schlieren photographs and hot-wire studies made by A. Demetriades 

also support these conclusions. Transition at 3 - 4 diameters in the 

P = 85 case is indicated, and at 4 - 5 diameters in the P = 60 psig case, 
o o 

In both cases it appears that the mechanism of transition commences in 

the 2 diameter region, corresponding to the position of the necked-down 

portion of the wake. (See Figure 15. ) The other possible explanation for 

the lack of predicted scale between the K Q turbulent P curves, is that 

the drag affects the scaling in a manner which cannot be determined by 

the present data. (The appearance of in the boundary condition of 

w 

the turbulent case was shown in Section III. 1. ) The present data i3 
inadequate in that momentum thickness could not be properly computed 
because of the effect of the tunnel perturbation. 

III. 4. 3. Transverse Spread of Concentration 

The similarity of concentration profile shapes may be checked by 
plotting the normalized concentrations K/K versus the normalized y 

IZiCtA 

coordinate y /y , where y„ is the mean of the oositive and negative 

absolute values and y is the same mean taken at the half-maximum 

m 

point of concentration. The curves are thus fitted at two points and the 
scatter pattern of the remaining points checked for similarity. The plots 
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of Figures 25 and 26 are the result of such a procedure. It can be seen 
that, except for the x/B = 3 case at P Q =10, the curves indicate a 
similarity as to shape of the spreading concentration profile. The 
exception noted in the x/D = 3 curve is to be expected, since it is too 
close to the origin of the similar region and an anomaly appears in the 
data of Figure 11. The profile points taken at P q = 85 after a slight 
shift of the flow, plot equally well with the profiles at 9 and 15 diameters 

before the flow shift. Having thus established at least reasonable shape 

-A 

similarity and the decay with x~ , one may plot K/K versus 



f - 



j/R. 



Y (x-k o )/; 



y d(x - x o ) 



and obtain one curve for all profiles. This is done for the P o = 85 case 
in Figure 27, with the appropriate virtual origin for the flow shift. In 
addition, a mean, J , and standard deviation, o" t are computed from 
the curve faired through the points, and the two parameters are used to 
compute a Gaussian distribution for comparison with the faired curve. 
The mean is just the width of the half-profile base and the standard 
deviation is given by 



I a*> 

= * < E A-ma*> 



for discreet values of the abscissa. The chi-squared test of statistical 
theory could be applied to determine the percentage confidence in the 
normality of the distribution. It was not done for this case because the 
tunnel perturbation effect would alter the meaning of the result. 
Reasonable certainty of agreement with the theory, however, is indicated. 
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In spite of the tunnel flow difficulties, it can be stated in summary 

that the flow is approximately similar for x/D = 9 throughout the Reynolds 

number range investigated. The related incompressible literature 

indicates that much greater downstream distance is required for the 

4 

flow to achieve complete similarity . One possible reason for this 
difference is that the proper transverse length scale for the hypersonic 
"inner wake" is actually much smaller than the body diameter, so that 
a distance of 9 - 10 body diameters corresponds to a distance of many 
inner wake momentum thicknesses. 

III. 4. 4. Non-Similar Regions of the Flow 

Very little quantitative evaluation of non- similar flow regions 

has been presented in the literature. This is understandable in view of 

the fact that these are regions of rapid adjustment in which nonlinear, 

fluctuating effects take place in the transition from one state of natural 

flow to another. In the wake of the hypersonic cylinder, such regions 

are to be found in the narrow portion of the wake and in the near vacuum 

at the rear of the model. A typical profile observed in the wedge shaped 

cavity between the model and the neck of the wake is that of Figure 28. 

The trapezoidal shape is evident until the concentration profile achieves 

a triangular distribution at the neck, after which the statistical effects 

of diffusion alter the shape toward that of the error curve. The base 

and point of the triangle are last to achieve the characteristic normal 

distribution in the similar flow downstream. This effect is also observed 

17 

in the incompressible coaxial jet profiles of Forstall and Shapiro . The 
trapezoidal region corresponds to their jet potential core. 
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III. 5. Suggestions for Improvement of the Experimental Method and 
Scope of the Wake Investigation 

III. 5. 1. Improvement of the Apparatus 

Several suggestions have previously been made concerning the 
apparatus, but a few additional points could be made, after reviewing 
the data obtained. 



III. 5. 1. a. Model Construction 

The data reveals that porosity of the model was adequate, but 
improved uniformity would possibly remove some of the unsymmetric 
behavior of the profiles. A few variations in the flow were noted which 
are not accounted for by the tunnel perturbation. It is possible that 
these resulted from model surface condition effects on the drag. The 
surface conditions are affected by model tool marks and other defects 
and the oil problem in the tunnel. These factors can be reduced to a 
minimum if the oil problem can be solved and the same model and model 
orientation is used for each test. The concentration measurements are 
very sensitive to drag effects on the wake pressure. 

III. 5. 1. b. Probe Installation 

The probe is restricted in the present installation to about 27 
diameters of travel, after moving the model rearward. The traversing 
mechanism should be modified to give about 50 - 60 diameters of travel 
downstream. It should be provided with a positive locking system to 
prevent shift during measurements, and a micrometer scale. It should 
be small enough in cross-sectional area to provide easy starting of the 
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flow. (This was a problem in the present investigation, with a probe 
mechanism of about 5/8 inch circular cross- section. ) 

III. 5. 1. c. Helium Flow Metering 

The metering of helium flow in the present investigation was 
satisfactory only after it was realized that model back pressure caused 
variances in the meter reading for the same mass flow, if the small 
meter was used. The large meter provided the necessary stability to 
differences in model porosity, but it was necessary to read the scale 
very accurately. The reading could be improved by the installation of 
a light and magnifying glass on the meter, with the provision of a means 
of sighting the level of the ball through a reticle or gunsight arrangement 
to insure consistency of eye level. With this arrangement the helium flow 
could be turned off and on at will and readjusted quickly, which would 
reduce the usage of helium during periods of apparatus adjustment and 
check-out. 



III. 5. 1. d. Sampling Apparatus and Procedure 

There are two major difficulties with the present sampling system 
that should be corrected before extensive testing is resumed: 

(1) The times required to clear the previous sample and bring 
in a new one are unnecessarily long. 

(2) The accuracy of the pressure correction is not satisfactory 
for very low sample densities. 

For the first problem, there are two possible solutions: The 
first has been previously mentioned as the pumping of each sample to 
barometric pressure. This method has attendant difficulties with regard 
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to possible contamination of the sample. In addition it requires a valve 
sequencing and pumping method which may become complex. It does, 
however, also eliminate the second problem. The second method is to 
follow the procedure utilized in the present investigation, measuring the 
pressure defect after a single compression, but to use larger compression 
ratio and provide automatic sequencing wherever possible. The automatic 
sequencing is desirable for precision and alleviation of the natural 
tendency to allow extra time for each step of the procedure. It could be 
accomplished by a suitable system of relays and solenoid-operated 
valves. The time would be even more reduced if a sensitive recorder 
is used in place of the potentiometer and the pen dropped for each sample 
point, in the relay sequencing. The proper time constants could be set 
in the relays with a few preliminary measurements at each reservoir 
pressure. The system and valves would have to be leak checked, of 
course, at frequent intervals. A rough diagram of a suitable system is 
included as an Appendix. 

The second problem, of measuring the pressure accurately, can 
be handled by reading total pressure intermittently with sample readings 
on the same recorder plot, if the sample is not compressed. Compression 
is desirable, however, at the very low tunnel densities, to prevent the 
pressure correction from becoming as large as the sample reading. If 
a transducer were installed on the sample volume the pressure could 
still be read intermittently on the recorder, and the total pressure at 
the probe could also be plotted as a third trace. 

The proper combination of the above ideas should provide equipment 
which would operate to take sample points automatically after each 
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movement of the probe. The Appendix is an example of one such possible 
combination. Sample points could be read in times of less than one 
minute if equipment of this type could be developed, instead of the 
several minutes average per point required in the tests of the present 
experiment. 

In addition to these problems, the water vapor contamination is 
of some additional concern. Either high tunnel make up air should be 
utilized and a reference sample gathered from the tunnel circuit, or 
low make up air with the single reference sample as taken in the present 
tests. 



III. 5. 1. e. Instrument Specifications 

The potentiometer of the present experiment was of Leeds and 
Northrup manufacture and was sensitive to four decimal places in 
millivolt units. It was well suited to the tests. Three decimal place 
accuracy is sufficient, but stability of zero is important. A recorder 
utilized in place of the potentiometer would need to meet the same 
specifications. 

The ammeter utilized was of the standard laboratory variety, 
accurate to about one-tenth milli-amp. An accuracy of 0. 01 milli-amp 
is desirable for maintaining constant current through the cell bridge. 

The conductivity cell itself was fairly accurate and stable, but 
it should be cleaned at intervals and checked frequently for possible 
drift of zero or potential bias of the setting made to read the initial 
calibration curve. Known samples are indispensible to the proper 
employment of the conductivity cell, and good mixing equipment and 
proceedu re should be employed for their preparation. The current and 












•> 
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voltage used in the cell bridge should be as high as possible for sensitivity 
but low enough to prevent fusing of the paired resistance wires. 

III. 5. 2. Suggestions for Additional Coverage of the 

Cylinder Wake 

Any new equipment designed for concentration studies in the wake 

should be checked against the results of Figure 15, which is the plot of 

points taken with the model in stable flow. Following check-out of the 

equipment, the profile data of this investigation should be repeated in 

the stable flow. The new data should include K points taken for 

max 

various tunnel pressures in the turbulent Reynolds number regions. 
Extensive total and static pressure coverage of the same P Q and x/D 
points should be made to determine exchange coefficients. Extensive 
coverage of the necked region of the wake should be obtained to determine 
how initial conditions might affect the later similar flow. Comparison 
of the normalized profiles with error distributions and modified 
formulas from the empirical theory should be made. Variation of the 
helium flow is desirable to see what effects may result in the diffusion 
and pressure profiles. Another diffusing gas should be employed to 
determine the effects in this regard. A heavy gas for instance might be 
useful in determining the relative diffusion of mass and momentum as 
compared to the light gas of the present experiment. An extensive program 
of wake coverage as outlined will provide a good background for studies of 
all types in wakes and boundary layers, where binary ga3 components 
may be detected by the thermal conductivity method. 
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IV. CONCLUSIONS 

The several important conclusions which may be derived from 
the present investigation are as follows: 

(1) The thermal conductivity method is an important and useful 
one for diffusion studies in the hypersonic flow field, where one binary 
constituent is to be detected in another. The method is adequate in 
determining mass diffusion properties of a flow with a light gas injected 
as tracer or as coolant for the model. 

(2) The problem of low densities in the samples taken from 
hypersonic flows can be controlled by suitable calibration. The procedure 
can be improved over that of the experiment with a more sensitive 
measurement of the sample pressures. 

(3) Other improvement of the model, probe, and instrumentation 
will increase the accuracy of the method employed in this case to within 
one or two per cent error for concentrations on the order of one-half 
per cent. 

(4) The method is adequate to determine turbulent processes of 
the flow and to detect transition, if the data is properly reduced and 
plotted to bring out the similarity properties. 

(5) The wake of the hypersonic cylinder is formed from the 
boundary layer gas which is retained initially in shear layers at very 
low pressure behind the model, in trapezoidal shaped profiles of mass. 

The gases are compressed into a very narrow wake at a distance of 
about two diameters downstream of the model, after which the lateral 
diffusion rapidly approaches the similarity regime. 

(6) The similar pattern is set up in either laminar or turbulent 
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flow at between 5-9 diameters behind the model. This result is in 
contrast to the much greater distances required in low speed flow. 

(7) In the laminar case, the theoretically predicted scaling of the 

maximum concentration inversely as the square root of the supply 
pressure was satisfied in the ratio of the P^ = 10 and P q = 35 psig 
readings. In the higher pressure cases (P q = 60, 85), did not 

scale with l/P as an approximate consideration indicates. Either the 
dependence on the pressure is more involved in the turbulent case, or 
the fully turbulent flow was not yet reached at the highest Reynolds 
number of the experiment. 

(8) Diffusion profiles were affected somewhat by a tunnel flow 
perturbation, but are still observed to be typically Gaussian, as theory 
predicts, a few diameters downstream of the narrowed portion of the 



wake. 
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TABLE I 

SUMMARY OF GAS PROPERTIES 



(1) 


(2) 


(3) 


(4) 


(5) 


Gas 


Molecular 

W T eight 


Ratio of 
Conductivity 
to Air at 

o°c 


Temperature 
Coefficient of 
Conduction 
0 - 100°C, 
per °C 


Per Cent 
in Dry 
Air at 
Sea Level 
By Mass 


Air 


28. 97 


1.00 


.00265 




Carbon Dioxide 


44. 01 


.605 


. 00527 


.05 


Helium 


4. 003 


5.97 


. 00256 


trace 


Hydrogen 


2. 016 


7. 15 


. 00265 


trace 


Nitrogen 


28. 02 


.996 


.00264 


75.80 


Oxygen 


32. 00 


1. 013 


. 00303 


23.22 


Argon. 


39. 99 


. 685 


. 00311 


.94 


Water Vapor 


IS. 02 


. 725 (46°C) 


. 00540 


0 - 6. 0 



Column References: 

<2) Handbook of Chemistry and Physics, 1st edition 

(3), (4) Daynes, Gas Analysis by Measurement of Thermal Conductivity, 
Table 1, Cambridge University Press, 1933 

Aeronautical Handbook, Pratt and Whitney Aircraft Company, 
1957 



( 5 ) 
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TABLE II 

CONCENTRATION TEST DATA SHEET 



Run No. Date 



Type of Run 




o5 



3 

w 

« 

H 



H 

W 



CO 



a 

2 

H 

U 

D 

Q 

W 

PS 

< 

H 

s 

a 

o 



3 

H 

a 

w 

o 

a 

o 

u 



i 
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TABLE IV 
ERROR ESTIMATE 



(1) 

Item 


(2) 

Cause 


(3) 

Amount 


(4) 

Range of 
Error 


(5) 

Possible 

Remedy 


1 . 

helium 

flow 

setting 


low sensitivity 
of flowmeter 


t 3.5% 
of emf 

(mv) 


all readings 


more sensitive 
meter 


2. 

zero 

shift 


unstable T. C. 
cell, ammeter, 
and 

potentiometer 


- .05 

mv 


all readings 


stable sensitive 
instruments; 
monitor current 
closely; keep 
cell clean 


3. 

pressure 

calibra- 

tion 

error 


low sensitivity 
of pressure 
scale 


t 3% 
to 0 % 
of mv 


K = 5%at 
Ap/p bs =1-0 

to >p/p bs < • 3 


use micro- 
manometer on 
pump scale; 
measure Ap to 
. 1 mm Hg 



4. Miscellaneous small effects having less than 1% aggregate error 
for average concentration readings of the tests: 



a. 


temperature difference 


e. 


vibration 


b. 


probe setting 


f. 


helium impurity (from lines) 


c. 


model differences 


g* 


tunnel conditions off by 
- . 2 psi and 2°F 


d. 


oil in flow 



5. Intermittent errors (large ones can usually be detected): 

a. water vapor in flow caused by excess make-up air with 
insufficient drying 

b. model installation or helium line leakage after metering 

c. large contimination in conductivity cell 

d. bad leaks in sampling system 

Total possible error for a sample of 11T helium concentration at pressure 
ratios A p/p bo < . 9 is approximately 6 % . 
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TABLE V 

TEST REYNOLDS NUMBERS 



T s Z75°r = constant 
0 


» P3ig 


Re/in x 10 ^ * 


Re^ x 10* * 


10 


5. 52 


1. 66 


35 


11. 18 


3. 35 


60 


16. 82 


5. 05 


85 


22. 48 


6. 75 



* Based on atmospheric pressure of 14. 4 psi for the laboratory 
elevation, and cylinder diameter, D = 0. 30 inch 



LIST OF PRESSURE DATA 
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This table is continued on the following page. 



TABLE VI -- CONTINUED 
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V 

W) 

a 

&o 

c 

£ 

o 



<D 

a 

o 

"d 

v 

d 

d 



o 

u 



J3 

d 



to 

3 

H 



TABLE VI -- CONTINUED 
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"SB 








© 




a 








a; 


0? 


u 


w 

3 


2 

QQ 


f—4 

0 


05 


cn 


<D 








Cli 






CD 


•w 



vu 

CQ 


<L> 


05 


CO 




0 ) 


a 




0 

4-1 


1— 4 

O 

a 


s 


a 

d 


(0 


4J 


II 


II 



0 . ft 
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TABLE VII 

SUMMARY OF CONCENTRATION DATA 



All Runs: M 


s 5. 8; T 

o 


3 275°F; xn„ = 0. 00305 t 
rle 


. 0001 lb/min 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Run 


Fig. 


Date 


p 0 . 


Type of 


x/ D 


Remarks 


No. 


No. 




psig 


Run 




‘ 


1 


10 


11/23/59 


85 


Profile 


3 




2 


10 


11/30/59 


85 


Profiles 


9,5 




3 


10, 13 


12/3/59 


85 


Profile 


15 


rh*. . 00279 Ib/min 

ine 


4 


11, 15 


12/22/59 


10 


Profile 


3 


5 


11, 15 


12/23/59 


10 


Profiles 


9 


xh *s . 00307 lb/xnin 












15 


iii-j. *=t .00281 Ib/min 
He ' 












24 


6 


14 


1/29/60 


85 


*“max 


3,5 
6. 5 


Flow shifted slightly 
since runs 1-5 












3,9 




7 


14 


2/2/60 


85 


K max P tS 


9,12,15 




8 


14 


2/4/60 


85 


K max 


28 






12 






■| Profile 


20, 36 


rh-j. *=. 00300 lb/min 
for G x / D - 20 


9 


14 


2/29/60 


35 


P ts 


9, 5, 6, 5, 3 
4, 2, 12, 20 






12, 13 






Profile 


15 




10 


14, 15 


3/18/60 


85 


K „ pt 
max 

K „ pts 
max 1 


4 


p.uns 10-15: 


11 


15 


3/21/60 


10 


3, 7 


Basic flow shifted 
again since runs 6-9. 


12 


15 


3/22/60 


10 


K pts 


5, 20, 27, 1 


Model position now 








35 


211w>A 


5, 15, 20 


2. 6” aft of position 


13 


15 


3/23/60 


35 


K pts 

max r 


5,27, 9,2 


for runs 1-9. 








60 


2, 5, 9, 20 





* m^j e values refer to continuity check. See Section III. 4. 



This table is continued on the following page. 
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TABLE VII- - CONTINUED 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Run 


Fig. 


Date 


P o- 

psig 


Type of 


x/D 


Remarks 


No. 


No. 




Run 






14 


15 


3/24/60 


60 


K ^ pta 
max 


5 3 7 
15,27 












Profile 


1 






14, 15 




85 


K a pts 
max * 


9, 5, 3, 2, 

1. 5, 18, 27 




15 


14, 15 


3/25/60 


35 


K „ pts 
max 


3,1 










85 




1.7 












SQiLmr/c of CALC/ 7 
LLr/'f v \iw<: k;w; Cunne. 



Steel nozzle blocks-? .Probe 
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DIAGRAM OF TEST SECTION AND EQUIPMENT, LEFT SIDE LEG I -FACING STREAM 




Gow-Moc Instrument Co 
Madison , N . J . 

Model 9285 with wiring 
mod. 91 93 ( Te - II ) 
Tungsten filaments 

Dimensions of cell: 
2x2x25/8" 



5 feed- through connectors wired 
in full Wheatstone bridge. Fiber- 
glass insulating sleeve, 12 " leads. 

For chromatogrophy to 300°C 



Sample 

aT 



EE 



Bf 



Ref. 

— E 



B+ + Pot- • B* 



Somple 

E3~rr 



I 

r 



Ref. 



=TZ3 



I B 



Top View 




FIG. 3 



THERMAL CONDUCTIVITY CELL 



7o 



Schemotic Zero contro1 




Zero control 
I or 2 Ohm pot. 



r 



Mongonin shunt 
selected for zero 
control of 40% of 
scole on either side B + 
of zero bolonce 




+ Pot.- 



Somple 




Wire wound pot. 
or helipot 
B- 



Current control 
20 Ohm pot. 



B + 

6 volts recommended 



Wiring diogrom 
for 4 filoment 
tungsten bridge 



FIG. 4 DIAGRAM OF WIRING FOR THERMAL CONDUCTIVITY 



BRIDGE 
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FIG. 5 MASS CONCENTRATION OF He IN N £ VS. BRIDGE 



DEFLECTIONS 
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FIG. 6 CALIBRATION OF CONDUCTIVITY CELL FOR PRESSURE 



DIFFERENCE BETWEEN SAMPLE AND REFERENCE 
CHAMBERS 
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Concentration , K , m.v. at arbitrary position in wake 
of cylinder injecting helium 



FIG. 7 



LINEARITY OF CONCENTRATION WITH A SMALL 
CHANGE IN FLOW OF HELIUM 



< 







FIG. 8 -- SCHLIEREN PHOTOGRAPH OF POROUS CIRCULAR 
CYLINDER, M = 5. 8 at P - 60 psig, T q = 275°F 
(extraneous pattern causeH by oil in flow; model 
shape affected by chips in glass) 




FIG 9 -- SCHLIEREN PHO' OGRAPH Oi PC' US GIF UL F 

CYLINDER, M - 5.8 at P = ps- ■ T ZT* 

o 



J1 



1.6 



1.4 



Concentration , K , % 



1.2 



Runs I - 3 '• M : 5.8 
P 0 =85psig, T 0 = 275°F 







FIG. 10 HELIUM CONCENTRATION PROFILES IN THE WAKE OF 

A POROUS CIRCULAR CYLINDER 



8 *. 




FIG. II HELIUM CONCENTRATION PROFILES IN THE WAKE 



OF A POROUS CIRCULAR CYLINDER 



n 







FIG. 12 HELIUM CONCENTRATION PROFILES IN THE WAKE 
OF A POROUS CIRCULAR CYLINDER 
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FIG. 13 COMPARISON OF PROFILES BEFORE AND AFTER 
SLIGHT TUNNEL FLOW SHIFT 
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FIG. 1 4 PLOT OF K m(JX VERSUS x/ Q BEFORE AND AFTER 
TUNNEL FLOW SHIFT 



r 
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FIG. 15 PLOT OF MAXIMUM CONCENTRATION VERSUS AXIAL DISTANCE 
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a. Comparison of flow for P 0 = 85 , x /D = 5 , model f wd., 
probe centered 





FIG. 16 FLOW SHIFTS WITH MODEL ATFWD. POSITION 
COMPARED WITH STABLE PLOT (C) 
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™He = 0 




FIG. 17 COMPARISON OF TOTAL PRESSURE TRAVERSES 
WITH MODEL IN STABLE AND SHIFTING FLOW 





b. P 0 = IO,x/D = 5, model fwd., probe at z=lin. 
off center 




c. Model at stable rear pos i t i o n , P 0 = 10, x/ D = 5 



FIG. 18 COMPARISON OF FLOWS FOR STABILITY AND 
MODEL POSITION 



2 



P/ 

P X P 

o 

xIO 2 



0 



0 




1. m He = 0 

2. m He = 0.003 lb. /min . 



I 



y/D 



-i.o o i.o 

a. With and without Helium flow , model fwd., probe 
centered 




b. With and without Helium flow, model fwd., probe 



at z = 1.0 




c. Duplicate trace for brass, porous models 
FIG. 19 FLOW COMPARISONS AT P 0 =85 AND 



x / D = 5 




FIG. 20 TOTAL PRESSURE AND MACH TRACES CORRESPONDING 
TO CONCENTRATION PROFILES, FIG. 10 
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FIG .21 



TOTAL PRESSURE AND MACH TRACES CORRESPONDING 
TO CONCENTRATION PROFILES , FIG. 1 1 
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FIG. 22 CONTINUITY PROFILE FOR HELIUM MASS PER 



SECOND IN THE WAKE OF THE POROUS CYLINDER 
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FIG. 23 PLOT TO DETERMINE VIRTUAL ORIGIN OF 
SIMILAR FLOW 




D 



FIG. 24 PLOT OF MAXIMUM CONCENTRATION VERSUS AXIAL 
DISTANCE WITH VIRTUAL ORIGIN CORRECTION 



1.0 



0.9 



0.8 



0.7 



r. 



96 



0.6 



k/k 



max. 



0.5 



0.4 



0.3 



0.2 



0. 



A V 






I 



6^ 



A 



A V 



/S 

A 



7 

Q 



Runs 1 , 3, 8 , 9 : 

M = 5.8, F^= 85 psig , T Q = 275°F 

m He = 0.0031 lb. /min. 

x/D : 
o 9 

& 15 Run 3 

15 Run 9 (after flow 
shift ) 

a 20 
? 36 






A 



£ 



gA 



3 A; 



□ L 



% 



7 
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FIG. 25 PLOT OF NORMALIZED CONCENTRATION VERSUS 



AVERAGE NORMALIZED VERTICAL DISTANCE 
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Runs 4,5: 

M = 5-8 , P Q = 10 psig , T q - 275° F 
m He = 0. 0031 I b. / m i n . 
x/D = 

o 3 ( non similar ) 
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FIG. 26 PLOT OF NORMALIZED CONCENTRATION VERSUS AVERAGE 
NORMALIZED VERTICAL DISTANCE 
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Gaussian OistriDution 
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where cr = 0.093 is the 
standard deviation of the 
faired curve with mean, 77 = 0 . 3 




M = 5.8 

P 0 = 85 ps ig 

T 0 = 275 °F 

m He = 0.0031 lb. /min. 

x/D : 

o 9 
A 15 

(Solid line faired 
through points) 



0 



0 . 



0.2 



0.3 



7 = v 4 dot^ 

FIG. 27 PLOT OF NORMALIZED CONCENTRATION TO SIMILAR 
COORDINATES VERSUS THE CORRESPONDING GAUSSIAN 
DISTRIBUTION 
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FIG. 28 TYPICAL PROFILE IN THE WAKE CAVITY NEAR 
THE MODEL , x/ D < 2 
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APPENDIX 

THERMAL CONDUCTIVITY APPARATUS AND PROCEDURE 
FOR HYPERSONIC WIND TUNNEL EMPLOYMENT 

Refer to Figure A-l 

I. Preliminary, with probe at desired distance behind model; 

1. Check all lines, volumes and valves for leakage. 

2 . Check concentration measurement on known sample. 

3. Start tunnel and stabilize. (Can be done along with 1 and 2 . ) 

4. Make zero concentration check with air in sample and 
reference cavities. Zero the concentration scale on the 
recorder. Turn helium on. 

5. Prepare to set pressure sensitivity by manually opening 
valve a, c, and d with m, b, and e closed. Measure total 
pressure on the micromanometer manually. Close d. Zero 
the pressure scale on the recorder. 

6. Set total pressure sensitivity in recorder ordinate scale; 
find tunnel centerline by probe movement vertically and set 
probe position potentiometer to read correctly on abscissa. 

II. Measurement of Total Pressure and Concentration: 

A. Automatic Operation, where sample is at pitot pressure 
greater than 1 cm Hg absolute; 

1. Set probe vertical distance, valves a, c, e, and f 

open, and b and d closed. Turn on "washing machine" 
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type relay system, which follows the following 
sequence: 

(a) Closes valve m and opens n to clear cell - 
relay holds for sufficient time, then opens 
m to allow next sample to wash through - 
relay again holds and then closes m. This 
operation leaves the ceil clean and at low 
vacuum for drawing of the tunnel sample. 

(b) Relay operates after sufficient time to achieve 
cell vacuum and closes n, then immediately 
opens m. (With l cm Hg or more differential, 
a tunnel sample is drawn into the cell in 
satisfactory time for automatic operation. ) 

After sample is drawn, relays switch to 
transducer circuit and drop recorder pen 
momentarily for total pressure, then switching 
to cell bridge circuit, drop pen for concentration. 

(c) Relay cycle is complete. Move probe and repeat 
for other points. This procedure obtains both 
total pressure and concentration. 

B. Manual Operation, where sample is at pitot pressure less 
than 1 cm Hg absolute: 

The pressure correction in this range is too great for 
direct concentration measurement of the sample at pitot 
pressure. The sample must be compressed. The hand- 
operated pump is adequate in that very few readings are 
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taken at these low pressures v/ith the present tunnel 
densities. The me rcury pump itself may be designed as a 
micromancmeter, in which case it may be used for the 
reading of cell pressures and will be more convenient for 
the manual operation than the transducer. (For a description 
of the hand operation, see Sections XL 7. 1. and III. 5. of the 
basic report. ) 

The measurement of total pressure along with compressed 
sample pressure requires additional steps. Usually, however, 
it will not be necessary to obtain total pressure at these low 
sample pressures, as they occur only in non-similar regions 
near the model. 

If possible, the cell bridge current should be turned off 



during each interval when the cell is completely evacuated. 
(See Section II. 7. 1. of the basic report. ) 
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® m.n : solenoid- operated valves 

Mercury pump 
(manual ) 



FIG. A— I : SCHEMATIC OF A SYSTEM TO MEASURE 

CONCENTRATION AND TOTAL PRESSURE AT LOW 
SAMPLE DENSITIES 





